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\ Abstract 

^""^ , We present in this paper a detailed note on the computation of Puiseux series solutions 

' of the Riccatti equation associated with a homogeneous linear ordinary differential equa- 

tion. This paper is a continuation of pQ which was on the complexity of solving arbitrary 
ordinary polynomial differential equations in terms of Puiseux series. 

Introduction 

IT) , 

{Sj \ Let K = Q(Ti, . . . , TJ) [77] be a finite extension of a finitely generated field over Q. The variables 

Ti,...,Ti are algebraically independent over Q and 77 is an algebraic element over the field 
Q(Ti, . . . , T[) with the minimal polynomial 4> 6 Z[T\, . . . , Ti][Z]. Let K be an algebraic closure 
of K and consider the two fields: 
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O- L = \J v& *K{{x»)) i C = U veW K{{x«)) 

which are the fields of fraction-power series of x over K (respectively K), i.e., the fields of 
Puiseux series of x with coefficients in K (respectively K). Each element ifi € L (respectively 
ifi € C) can be represented in the form tp = YlieQ c i x ^ Ci £ K (respectively q € K). The 
order of ip is defined by ord(ip) := min{i £ Q,q 7^ 0}. The fields L and C are differential 
fields with the differentiation 

Let S(y) = be a homogeneous linear ordinary differential equation which is writen in the 
form 

S(y) = f n y {n) + ■■■ + fiy' + f y 

where fi G K[x] for all < i < n and f n ^ (we say that the order of S(y) = is n). 
Let yo,...,y n be new variables algebraically independent over K{x). We will associate to 
S{y) = a non- linear differential polynomial R € K [x] [yo, . . . , y n ] such that y is a solution of 
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S(y) = if and only if ^ is a solution of R(y) = where the last equation is the ordinary 

differential equation R(y, . . . , = 0. We consider the change of variable z = i.e., 
y' = zy, we compute the successive derivatives of y and we make them in the equation 
S(y) = to get a non-linear differential equation R(z) = which satisfies the above property. 
R is called the Riccatti differential polynomial associated with S(y) = 0. We will describe all 
the fundamental solutions (see e.g. [201 [T3]) of the differential equation R(y) = in £ by a 
differential version of the Newton polygon process. There is another way to formulate R: let 
{fi)i>Q be the sequence of the following differential polynomials 

r = 1, r 1 =y ,..., r i+1 = y rj + Dr u Mi > 1, 

where Dyi = y.- L+ i for any < i < n — 1. We remark that for all i > 1, n £ ^[yo 5 • • • >yi-i] 
has total degree equal to i w.r.t. yo> • • • > yj_i and the only term of ri of degree i is y Q . 

Lemma 0.1 The non-linear differential polynomial 

R = frJn H V fin + foro G #M [yo, • • • , Vn] 

is the Riccatti differential polynomial associated with S(y) = 0. 

1 Newton polygon of R 

Let R be the Riccatti differential polynomial associated with S(y) = as in Lemma lO.ll We 
will describe the Newton polygon M(R) of R in the neighborhood of x = +oo which is defined 
explicitly in [lj. For every < i < n, mark the point (deg(/i),i) in the plane M 2 . Let be 
the convex hull of these points with the point (— oo, 0). 

Lemma 1.1 The Newton polygon of R in the neighborhood of x = +oo isN ' , i.e., Af(R) = N . 

Proof. For all < i < n, deg yo = i and the only term of of degree i is yg, 
then lc(fi)x deg ^yQ is a term of R and Af C J\f(R). For any other term of fri in the form 
bxiyQ ■ ■ ■ y£~x , where b € K, j < deg(f) and oto + ■ ■ ■ + Oj-i < i, the corresponding point 
(j — a\ — ■ ■ ■ — (i — l)aj_i, «o + • • • + oti-i) is in the interior of M. Thus M(R) C M. □ 

Lemma 1.2 For any edge e of N{R), the characteristic polynomial of R associated with e is 
a non-zero polynomial. For any vertex p ofAf(R), the indicial polynomial of R associated with 
p is a non-zero constant. Moreover, if the ordinate of p is io, then hrjt tP \(fj,) = lc(fi ) ^ 0. 

Proof. By Lemma [TTT1 each edge e S E(R) joints two vertices (deg(/j 1 ), i\) and (deg(/j 2 ), i-i) 
oiM(R). Moreover, the set N(R,a(e),b(e)) contains these two points. Then 

+ H {R>e) (C) = /c(/ n )C il + lc(f l2 )C^ + t. 

where t is a sum of terms of degree different from i\ and ii- For any vertex p € V(R) of 
ordinate io, /c(/j )x deg ^ i o)yQ° is the only term of R whose corresponding point p. Then 

h(R, P )(v) = Hf k) ) ^o.n 

Corollary 1.3 For any edge e G E(R), the set Am e \ is a finite set. For any vertex p £ V(R), 
we have Am p \ = 0. 
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2 Derivatives of the Riccatti equation 

For each i > and k > 0, the fe-th derivative of rj is the differential polynomial defined by 



(0) ._ (l) ._ / ._ (fc+i) ._ ( (fc) y _ g rj 



Lemma 2.1 For a// i > 1, we have r[ = ir^—i. Thus for all k > 0, r- = (*)jfe7"i jfe, where 

(i)o := 1 and := i(* — 1) ■ • • (i — A; + f). 

Proof. We prove the first item by induction on i. For i = 1, we have r^ = 1 = l.ro- Suppose 
that this property holds for a certain i and prove it for i + 1. Namely, 

= 2/or- + ^ + £>r- 

= iyo^i-i + n + D(in-i) 

= i(yon-i + + rj 

= in + rj 

= {i + i)n. 

The second item is just a result of the first one (by induction on A;). □ 

Definition 2.2 Let R be the Riccatti differential polynomial associated with S(y) = 0. For 
each k > 0, the k-th derivative of R is defined by 

* {k) ■= S = E M k) - 

° y 0<i<n 

Lemma 2.3 For all k > 0, we raaue 

0<i<n-fc 

Proof. For all i < k, we have rj fc ^ = because deg OT (rj) = i. Then by Lemma I2TT1 we get 

R {k) = E M k) 

k<i<n 

= E fi( i h r i-k 
k<i<n 

= E & + k )kfj+kTj 
0<j<n-k 



with the change j = i — k. □ 

Corollary 2.4 T/ie &;-£a derivative of R is the Riccatti differential polynomial of the following 
linear ordinary differential equation of order n — k 



S {k) (y):= E (*' + *0*AW 



0<i<n-A: 

Proof. By Lemmas 10.11 and 12.31 □ 
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3 Newton polygon of the derivatives of R 

Let < k < n and RV*' be the fc-th derivative of R. In this subsection, we will describe the 
Newton polygon of R^ k \ Recall that R^ is the A;-th partial derivative of R w.r.t. yo> then 
by the section 2 of pQ , the Newton polygon of R^ is the translation of that of R defined by 
the point (0, -k), i.e., Af(R^) = M(R) + {(0, -k)}. The vertices of A/"(i? (/c) ) are among the 
points (deg(/j + /c), i) for < i < n — k by Lemma [231 Then for each edge of Af(R^) there 
are two possibilities: the first one is that is parallel to a certain edge e of Af(R), i.e., is 
the translation of e by the point {(0, — k)}. The second possibility is that the upper vertex 
of ek is the translation of the upper vertex of a certain edge e of Af(R) and the lower vertex 
of e/c is the translation of a certain point (deg(/j ), io) of Af(R) which does not belong to e. 
In both possibilities, we say that the edge e is associated with the edge e^. 

Lemma 3.1 Let eu £ E(R^) be parallel to an edge e € E(R). Then the characteristic 
polynomial of R^ associated with e^ is the k-th derivative of that of R associated with e, i.e., 

H (RW,e k )(. C ) = H (R,e)( C )- 

Proof. The edges e^ and e have the same inclination /x e = [i eh and N(R( k \a(ek),b(ek)) = 
N(R, a(e), 6(e)) + {(0,-fc)}. Then 

H (R^,e k )(C) = Yl (i + k) k lc{h+k)C* 

(deg(f l+k ),i)eN(RW,a(e k )Mek)) 
(deg(/,)J)e7V(i?„a(e),fe(e)) 

Corollary 3.2 For any edge e k £ E(R^), the set A, R (k) 6fc ) is a finite set, i.e., H^ R ( k ) efc )(C) 
is a non-zero polynomial. For any vertex pk € V{R^), we have A^ R ( k ) Pk \ = 0. 

Proof. By Corollaries ELl and □ 

4 Newton polygon of evaluations of R 

Let R be the Riccatti differential polynomial associated with S(y) = 0. Let < c £ K , /i £ Q 
and R\(y) = R(y + cx^). We will describe the Newton polygon of Ri for different values of c 
and [i. 

Lemma 4.1 R\ is the Riccatti differential polynomial of the following linear ordinary differ- 
ential equation of order less or equal than n 

Si(y) ■= £ ±R®(ca?)y®. 

0<i<n 
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Proof. It is equivalent to prove the following analogy of Taylor formula: 

Ri= J2 \R {i) {^)n 

0<i<n 

which is proved in Lemma 2.1 of [13J . □ 

Then the vertices of Af(R\) are among the points (deg(R^ (cx^),i) for < i < n. 
Thus the Newton polygon of R± is given by (Lemma 2.2 of |13j): 

Lemma 4.2 If fi is the inclination of an edge e of M{R), then the edges of M{R\) situated 
above e are the same as inM(R). Moreover, if c is a root of H( Re ^ of multiplicity m > 1 then 
Af(R\) contains an edge e\ parallel to e originating from the same upper vertex as e where 
the ordinate of the lower vertex of e\ equals to m. If m = degH( R ^, then J\f(R\) contains 
an edge with inclination less than jj, originating from the same upper vertex as e. 

Remark 4.3 // we evaluate R on cx^ we get 

R{cx») =^/ iX (cV" + t), 

0<i<n 

where t is a sum of terms of degree strictly less than Then 

lc{R(cx»)) = YsHhY = E = H (RM C )> 

i&B (deg(/i),i)ee 

where 

B := {0<i<n; deg(/j) +ifi= max (deg(/ j ) + j>; f j ^ 0)} 

0<j<n 

= {0 < i < n; (deg(/i),i) € e and fi / 0}. 

Lemma 4.4 Let /i be the inclination of an edge e of Af(R) and c be a root of Hm e ) °f 
multiplicity m > 1. Then 

H(R ltei )(C) = H( Re )(C + c) 

where e\ is the edge of M{R\) given by Lemma In addition, if e' is an edge of M{R\) 
situated above e (which is also an edge of N{R) by Lemma \4Jfy then H( Rl ^(C) = H( R ^(C). 

Proof. We have 

H {R , e) (C + c) = £ ^ e) (c)C* 

m<k<n 
m<k<n 

= ^HR {k) {cx»))C k 

m<k<n 
= ^(iii,ei)(C) 

where the first equality is just the Taylor formula taking into account that c is a root of Hm e \ 
of multiplicity m > 1. The second equality holds by Lemma f3.1[ The third one by Remark l4.31 
The fourth one by Lemma 14.11 and by the definition of the characteristic polynomial. □ 
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5 Application of Newton-Puiseux algorithm to R 



We apply the Newton-Puiseux algorithm described in pQ to the Riccatti differential polyno- 
mial R associated with the linear ordinary differential equation S{y) = 0. This algorithm 
constructs a tree T = T(R) with a root tq. For each node r of T, it computes a finite field 
K T = K[9 T ], elements c T G K T /i T £ QU {— oo,+oo} and a differential polynomial i? r as 
above. Let IA be the set of all the vertices r of T such that either deg(r) = +00 and for the 
ancestor t\ of r it holds deg(ri) < +00 or deg("r) < +00 and r is a leaf of T. There is a 
bijective correspondance between IA and the set of the solutions of R{y) = in the differential 
field C. The following lemma is a differential version of Lemma 2.1 of f4j which separates any 
two different solutions in C of the Riccatti equation R(y) = 0. 

Lemma 5.1 Let ipi,ip2 € C be two different solutions of the differential Riccatti equation 
R(y) = 0. Then there exist an integer 7 = 712, 1 < 7 < n, elements ^1,^2 £ K, £1 7^ £2 
a number € Q such that 

ord(R^(^i)-^ 12 ) < Mia, /«" = 1,2. 

Proof. By the above bijection, there are two elements u\ and U2 of IA which correspond respec- 
tively to ipi and tp^- Let i$ = max{i > 0;rj(ui) = Tiivq)}. Denote by r := Ti {u\) = Ti (1*2) 
and ti := Tj 0+ i(ni), T2 := Tj 0+ i(«2). We have n / T2 and e := max(/U n , /i T2 ) is the 
inclination of a certain edge e of M(R T ). There are three possibilities for e: 

- If /u T2 < |U ri then e = /i ri = /i e . We have c T1 is a root of H( Rr ^ of multiplicity mi > 1 
and R Tl = R T {y + c^x^i). Then by Lemma |4~21 there is an edge e± of Af(R Tl ) parallel to e 
(so its inclination is e = [i T1 ) originating from the same upper vertex as e where the ordinate 
of the lower vertex of e\ equals to m\. In addition, e is also an edge of Af(R T2 ) and by 
Lemma WM we have H^ Rt >e )(C) = H^ RtS ^(C) and 

^ 1 ,e 1 )(C)=^,e)(C + C Tl ). (1) 

- If fi n < fi T2 then e = /i T2 . Then by Lemma POI there is an edge e2 oiM(R T2 ) parallel to e 
originating from the same upper vertex as e. By Lemma f4.44 we have e ) (C) = H( Rt ^ (C) 
and 

H {RT2fi2) (C) = H (RTte) {C + c T2 ). (2) 

- If /i Tl = fi T2 = e then c Tl and c T2 are two dictinct roots of the same polynomial 
H(R Tt e)(C). Then equalities of type ([I]) and ((2]) hold. 

Set 7 = deg<7(-£f(fl Tie )) — 1 < deg yo Vn (R) — 1 < n — 1 < n and 7 > 1 because that 
the polynomial H( Rr ^(C) has at least two distinct roots c n and c T2 . Moreover, we have 
ord(ipi — y Ti ) < e for i = 1,2. Let £j € K Tl and £ 2 £ -^r 2 be the coefficients of C 7 in the 
expansion of H^ Rt ^^(C) and H^ Rt ^{C) respectively. There is a point (//i2,7) on the edge 
e which corresponds to the term of H^ Rt ^(C) of degree 7. We know by Lemma |4. II that 

R(y + ipi)= ^ RU) ^i) r 3 for » = 1, 2. 

0<i<n ^' 

Then ord{R^\ipi) - 7!$ jX^ 12 ) < /J12 for i = 1,2. This proves the lemma by taking £j = 7!^ 
for i = 1,2. □ 
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Let {^>i, . . . ,^/ n } be a fundamental system of solutions of the linear differential equa- 
tion S(y) = (see e.g. [201 Hffl [33]) and tpi, . . . ,i(j n be their logarithmic derivatives 
respectively, i.e., Vi = ••• = ^n/*n- Then = for all 1 < i < n. 

Definition 5.2 Lei -0 ^ e ara element of the field C We denote by span r {ip) the r- differential 
span ofip, i.e., span r {ip) is the Z-module generated by ri(ip), r2(tp), ■ ■ ■ ■ 

Lemma 5.3 Let ip E C be a solution of a Riccatti equation i?2(y) = where R2 € 
Z[yo, . . . , y n ] of degree n. Then span r (ip) is the Z-module generated by ri(ip), . . . , r n -i(ip). 

Proof. Write R2 in the form R2 = r n + a n _ir n _i + • • • + a±ri + o.q where ai € Z for all 
< i < n. Then 

r n+1 (ip) = tpr n (ip) + Dr n (ip) 

= ai(ipri(i/)) + Dn(ip)) 

0<i<n 
0<i<n 
0<i<n 

for suitable fa G Z using the fact that r n (V>) = X^o<j<n a i r i(' l P)- a 

Consider a Z-module M := span r (ipi, . . . ,ip n ), i.e., M is the Z-module generated by 
ri(^i),r2(ipi), ■ ■ ■ for all 1 < i < n. We define now what we call a r-cyclic vector for M (this 
definition is similar to that of the cyclic vectors in [161 \§\ 119]). 

Definition 5.4 An element m G M is called a r-cyclic vector for M if M = span r {m). 

The following theorem is called a r-cyclic vector theorem. It is similar to the cyclic vector 
theorem of [El H EJ [19] . 

Theorem 5.5 Let M be the "L-module defined as above. There is a r-cyclic vector m for M . 

Corollary 5.6 Let m £ M be a r-cyclic vector for M . Then for any 1 < i < n, there exists 
a Riccatti differential polynomial Ri S Z[yo, • • • ,Vn] such that ipi = Ri(m). 

Lemma 5.7 For each element m G M , one can compute a Riccatti differential polynomial 
R m £ if [x] [7/0, ■ ■ ■ ,Un] such that R m (m) = 0. In addition, there is a positive integer s such 
that the order of R m (y) = and the degree of R m w.r.t. y$, . . . ,y n are < n s . 

Proof. Each element m 6 M has the form m = aitpi + • • • + a n ip n where a\, . . . , a n £ Q. 
Then 

m = 111 n ' 

is the logarithmic derivative of • • • Or Lemma 3.8 (a) of [T7] (see also [151 IS]) 

proves that one can construct a linear differential equation S m {y) = 0, denoted by 
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S® ai+ '" +an (y) = of order < n ai+ - +a " such that ty" 1 ■ ■ ■ is one of its solutions. The 
equation +a " (y) = is called the {a\ + • • • + a n )-th symmetric power of the linear dif- 

ferential equation S(y) = 0. In order to compute the equation S m (y) = associated with the 

linear combination m = aiV'i H h ct n ip n € M, we take the change of variable z = y ai ^ ha ™ 

where y is a solution of S(y) = and we compute the successive derivatives of z until we 
get a linear dependent family over K. The relation between these successive derivatives gives 
us the linear differential equation S m (z) = 0. Let R m be the Riccatti differential polynomial 
associated with S m (y) = 0, then m is a solution of the equation R m {y) = 0. □ 



Remark 5.8 For any 1 < i < n, we can take Rj )i = R where R^ i is defined in Lemma 5. 7 
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